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INTRODUCTION 

The STARS programs perform a significant amount of calculations prior to the 
complete processing of all data. Thus data errors or inconsistencies usually 
cause ‘a run to terminate after the first data error is encountered.- To 
overcome the resulting problem of requiring several computer runs to debug 
data, special data debugging packages called SAT-1S, SAT-1B, and SAT-lV 
were cseated for the STARS statics, stability, and vibrations programs 
respectively." ' ’ ' 

In the analysis of shells of revolution subjected to unsymmetric loads, the 
loading is often available as point intensities, whereas the analysis programs 
require Fourier series input. The SAT-2S program provides the necessary 
conversions. This program could have been incorporated into STARS but. it 
was considered that the user should have additional freedom in specifying 
the series truncation points, and incremental summing options, all based 
upon his knowledge of the accuracy with which the loads were originally 
obtained. 

The programs were written exclusively in FORTRAN IV for the IBM 370-165 
computer, and, then converted to the UNIVAC 1108. The core utilization of the 
programs is: 

SAT -IS : 35,700 words of core 

SAT-2S: 26,700 words of core 

/ - 1 1 • * 

SAT-1B : 31,100 words of core 

c SAT- IV : 31,100 words of core 
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SECTION 1 

SATELLITE -IS PROGRAM 

> *j ' K ’t - ( f- — 


1.1 INPUT : The sole purpose of the SAT-l( ) series of programs is data , 

debugging of the appropriate STARS program input data decks. Therefore 
..the input data is identical to that described for the STARS programs, 
(Refs. 1,2), with one exception to be discussed below. For this ... 
purpose the dash separator cards, which are actually required only for 
data separation in the SATELLITE programs, are accepted also byivthe 
STARS programs. , . 

1 The one' data change involves the first card of the data deck, which in 
the STARS programs is an arbitrary alphanumeric card. For the SAT-l( ) 
series of programs this card is 


; I. Title Card Column 

A. i* Alphameric title (submission description) 1-60 
, B. Scale Extent 61-70 

This number will be used as the scale 
extent in the graphics output of the 
SATELLITE program. For example if the 
input number is 3000. , the scale for 
all the diagrams will be set (square 1 

scale is used) so that 3,000.0 fits •: 
properly on 8 in. paper. A "default"! ' 

option to this input is also available. 

If no scale extent is input, the program 
will size each region topology diagram 
on a square scale to best fit that 
region, resulting in different scales 
for different regions. 


Forma 
15A4 
F10 . 0 




All other input remains unchanged from the appropriate STARS program. 
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1.2 OUTPUT ; The SAT-l( ) program functions basically are the following: 
...A,. . Check all STARS input data for' consistency. , . t 

B. ..(Check all STARS input data for. possibilities of causing . "divide 

by zero" errors. • • • . . • 

C. Check shell idealizations by providing exact plots of 'input points 
and/or shapes, and topology, thus allowing the user to locate 
gaps or other types of errors. 

4 . 

4 

- V' D. ' Check- spellings on alphameric clues, and where ^possible, format 

errors. In the latter case, no attempt has been made to override 
systems automatic terminations due to format inconsistencies. 
However, the use of the dash separator cards often allows further 
checking to proceed. The dash separator cards also overcome user 
errors in setting input table lengths. ■■ 

Many of the errors encountered in the STABS data deck will not 
affect the idealization plotting capability of the SAT-l( ) 
programs. In cases where plotting is affected, as much plotting 
as possible will be accomplished before termination. Additional 
error, messages will be provided directly on the charts. 


1.3 EXAMPLES, FLOW CHART, LISTING ; Two sample test problems are executed 
by the SAT-1S program. The first problem involves explicit input 
errors in a multi -harmonic submission, while the second test problem 
error is in idealization. The first problem errors are: 

f 1. In the first segment, first region geometry table the following 
items are all zero: iso-grid spacing, the inner and outer sheet 
thicknesses of the basic sandwich shell wall, the core thickness 
of the basic sandwich shell wall. 

2. The wall geometry clues are in error. There cannot be a ring in 
the meridional direction (first clue), and the second clue is 
misspelled. 
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3 . 

■'4.' 

- 3 ; 5 -.' 

6 . 

I.‘‘ . 

i 

T. 

8 . 

-'9. 

10 . 


In the second segment, first region geometry table the following 
'items are all zero: waffle grid spacing, the inner ;and' outer 

sheet thicknesses of the basic sandwich shell wall,; the core 
thickness of the basic sandwich shell wall. . 

The second wall geometry clue for segment two,: region :one, is 
misspelled. • • nvv 

The joints for the intra-region link in region one are-reversed. 

In the first segment, second region geometry table the following 
items' are all zero: K ll' K 1 2’ K 2V K 22’ K 33’ D ll’ D I 2’ D 21’ D 22’ 

ar - 4 v ' . . 'IZ'Z 

The above error (item 6) is repeated' . in , the second segment of 
the second region. 

The. structure kinematic link ; cards are. reversed. 

Because of the kinematic linking, the degrees of freedom of 
joints 1 and 2 must be "zeroed out" . This run is obviously 
a fictitious test run since all joints cannbt be "zeroed out", 
and thus joint 3 must have some degrees of freedom remaining. 

The error count is printed as 63. 
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FOR INFORMATION Call V. SVALBONaS 



NOTE - GRAPHICS output WILL. Be produced 




the iso-grid spacing is zero 



. BE-GlQN- NUha ER L_±lliE«MENT:.NU«aM l SEGMENT GEOMETRY TABLE IlEM 



region number i segment number i segment geometry table item 



REG.IQ.N..JLUBBER .1 SEGMENT NUMBER 1 SEGMENT GEOMETRY table item 
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on number i segm ent number 2 segment geometry, table iteh 



GION number i SEGME N T NUMbe R - 2 SEGMENT . GEOMETRY J able . ITEM. 



REGION NUMBER 1 SEGMENT NUMBER 



region Dumber 2 segment number' i segment' geometry table item 




region number * segment number i segment geometry table item 



THt._0ai._5T i femes' 5 parameter is zero*. 



regTo'n* number •••‘jj— tt- segment number i segment geometry table item 



region number 2 segment number 1 segment geometry tabie item 



the eh stiffness parameter is zero 



REfiJ.eN.-NU.NBER .SEGMENT number 2 SEGMENT GEOMETRY Table item 



region number 2 segment. number 2 segment geometry table item 





..REGION NUM«£ft Z SEGMENT NUHBER 2 SEGMENT GEOhCTRY TABLE ITEM 




20000C0*DH ^cqooqo+o* 







.0000000 .785*4000*00 >1570800*01 .2356200*0) .3141600*01 

• 3V27 000*0 1 ^47+2400*01 "5497800*01 




im M' ‘fyoCHi pflj* i 

REGION NUMBER 1 2 SEGMENTS 1 LINKS 
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JOJ AlP H?00**> PAJ f ? 
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The second problem has no detectable error printout. However, inspection of 
the graphical display of the idealization shows an unintentional wrinkle put 
in by the point input geometries. This means that the z, r tables were 
acceptable to the program, but not correct for the representation of a 
smooth, toroidal-like shape. This type of input error would not cause the 
STARS program to interrupt execution, but would produce an analysis of the 
wrong structure. 









for information call v. svalsonas 



note * graphics output will be produced 
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The flow charts and listings of the SAT -IS program follow. The functions of 
the MAIN, GRAPH, ZPIGT and BOZO routines are explicit from the flow charts 
and need no further elaboration. The routines GEOMET, PLINE and PLICO are 
discussed in References 1, 2. 
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NERR0R = 2 
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RE AD (5*110) JN,DLP,ANGL,DUM 
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SECTION 2 


SATELLITE -2S PROGRAM 


2.1 INPUT: The SAT-2S program takes a point input and obtains .a Fourier 

series representation for the points. Thus such loading as that 
caused by maneuvers (aerodynamic loading) which is usually known as 
point intensities, can be converted to STARS-2S Fourier input. The 
decision to create a separate SATELLITE program for this function, 
rather than a subroutine in STARS-2S was based on the following 
considerations: 

(i) Only the user knows the accuracy with which the loads were 
calculated. Thus he alone is best able to decide the accuracy 
of the Fourier series representation required . 

(ii) By checking the graphical output in SAT-2S the user can 
determine whether he has input enough points to accurately 
represent his load distribution in Fourier series form. 

The user is also able to see the Fourier series convergence 
trends and make better decisions on truncation. 

(iii) The user is able to specify incremental summing of the 
Fourier series, thus eliminating STARS program runs with 
harmonics which are infinitesimal (but not zero) compared 
to the dominant harmonics. 


The order of input for the SAT-2S is the following. 


I. Program Control Card 

A. Number of input points 

The number of input points must be 
between 20 and 720. As a general 
rule for rapidly varying distributions 
("almost" point loads) a dense mesh of 
input points should be provided. A good 
representation is about a 12 to 15 point 
density per harmonic full wave. 


Column 


1-3 


Format 

13 
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B. Type of series identification 

0 = full Fourier series .required 

1 = cosine half series required 

2 = sine half. series required 

G. First truncation harmonic 7-9 13 

This is the first point to which the 
series will be summed. 

D. .Second truncation harmonic 10-12 13 

This is the second point to which the 

series will be summed. 

E. Total number of harmonics required ^ 200 13-15 13 

This is the total number of series terms 

which will be calculated. The final 
series will be summed to this number of 
terms . 

F. Graphics clue l6-l8 13 

0 = no series plots 

1 = include series plots 

G. Series summing increment or multiple (INC) 19-21 13 

INC = 0, 1, -1: sums will involve all terms 

INC > 0: sums will consist of the following 
terms: zeroth + (INC) + (2INC) + 

(3INC) + 

INC < 0: sums will consist of the following 
terms: zeroth + first + (l+INC) + 

( 1+2 INC ) + (1+3INC) + 


Column Format 

k-6 13 
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Column 

Format 

II. Load Point Intensity Cards 

A. Point intensities (between 20 and 720 ) 

1-70 

5E14.7 

Use as many cards as necessary. The 
point intensities are assumed to be 



equally spaced, by the program, between 
0° and 360° . The first input intensity 
should be that at 0^. & nd the last should 
be that at 360° . The program will set 
the 0° intensity value equal to that at 
360° . Tbe intensities are signed (+ or -) 

values. 



2.2 OUTPUT : The SAT-2S program output is the following: 

A. Title Page - The title of the program is printed along with 
pertinent information from the Program Control Card. 

B. Harmonic Amplitudes - The Fourier series sine and cosine 
coefficients are printed next. 

C. Fourier sums - The input points and the Fourier series 
representations are printed next. The increment value (if any 
used) is printed among the headings. The Fourier series will 
be summed at the first two input locations, and finally the 
total sum will be printed. 

D. Graphics - If requested the graphical output will display the arrays 
of printed output (item C) on three charts - one for each truncation 
point. The input points are displayed as a solid line, and the 
Fourier series approximations as x's. The Fourier series 
approximations are charted at every degree (360 points) regardless 
of the input point number. In this way unexpected irregularities 

in the series between input points can be noted. 
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2 . 3 EXAMPLES, FLOW CHART, LISTING : A sample test problem is executed 

by the SAT-2S program. The loading considered is made up of two 



The program is required to find a 90 harmonic Fourier series 
representation of the load, and also to print a 30 and 60 
harmonic estimate. 

From the output it can be seen that all even terms are insignificant 
and INC could have been set equal to -2. The even terms should not 
be used in a STARS run. The graphical output exemplifies the fact 
that about 12 to 15 points per harmonic are necessary for an accurate 
calculation of the harmonic amplitude. The improvement in convergence 
is easily seen in increasing the number of series terms from 30 to 60. 
With a 720 point input, the 60th harmonic has a 720/60 = 12 point 
representation. Increasing the number of series terms to 90 improves 
convergence in the area of the spikes, but causes more eratic 
behavior around zero. The 90th harmonic has a 720/60 = 8 point 
representation. Thus the harmonic amplitudes from 6l to 90 are not 
calculated with appropriate accuracy. 
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The flow chart and listing for the SAT-2S program follow. There is only one 
routine, MAIN, aside from system graphics routines. 
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MA'P ‘,15 A, B 
LIB SYS**MSFC$ 
IN MAIN 
XOT B 



SECTION 3 


SATELLITE -IB PROGRAM ! 

The discussions in SECTION 1 for items 1.1, 1.2 are equally valid for all 
SAT-1 ( ) programs. 

EXAMPLES, PLOW CHART, LISTING : Two sample test problems are executed 

by the SAT-1B program. Both problems involve explicit input errors. 

The first problem errors are as follows: 

1. The ring specification number on the first region card exceeds 
28, the maximum number of rings allowable. 

2. The number of points input on the first segment exceeds l4, the 
maximum number allowable. Since the program does not know if 
the input number 15 is in error, or if the user actually input 
an extra point, it scans for the next separator card before 
restarting. After printing an appropriate message the cards 

up to the separator are not checked, but printed in card image 
form. 

3. The clue for the second segment should be either l^A or l4B. 

k. After finishing region one, the program scans the ring cards 
assumming the input number, 29, as being correct. It finds 
that on all 29 rings, the ring centroid radius is input as zero. 

5. The second region, second segment has the errors noted in items 
3 and 2 above. 

6. The structure ring specification number has also exceeded 28, 
and again when the ring cards are checked a zero ring centroidal 
radius is found on all rings. 

7. The error count is printed as 6b. 
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The second problem errors are as follows: 

1. The recommended sizing for an annular plate segment has been violated 

in the first' region, first segment ... 

2. The thermal clue THST is inconsistent with the first entry, "1", on 

the load clue card for the dead load. .,, r . .. . 

3. The first entry on the load clue card for the live load actually 
causes three errors: 

a) Both dead and live loads cannot contain temperature. 

b) The entry, "5", is inconsistent with the thermal clue THST. 

c) The entry, "5", is invalid, and causes the program to scan 
further cards for the dash -separator card. 

h. In the geometry table for the second segment, first region the following 
items are all zero: waffle grid spacing, core thickness of the basic 
sandwich wall, the , sheet thickness of the basic sandwich wall. 

5. The same thermal clue versus load clue errors are present in this 
segment as those discussed in items 2 and 3 above. The difference 
here is that the first entry on the live load clue card is "4", thus 
allowing a check of the load cards. 

6. A check of the thermal loads against the appropriate material property 

I 

table reveals that the table has an insufficient thermal range. 

7. In the first segment, second region the following geometry table items 
are all zero: ring spacing, stringer spacing, the sheet thickness of 
the basic shell wall. 

8. In the first segment, third region the following geometry table items 
are all zero: waffle grid spacing, the outside and inside sheet thicknes 
of the basic sandwich shell wall, the core thickness of the basic 
sandwich shell wall. 
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9. The structure kinematic link cards are in reverse order. 

10. Joint 1 must be zeroed out in all 1 the boundary condition sets, due 
to the kinematic link connection. 

11. The error count is printed as 33- 
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' : I 

! ! I 
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RADIUS OF REVOLUTION 
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The flow charts and listing of the SAT-1B program follow. As can be seen 
from the MAIN routine flow diagram the only major difference between the 
SAT-lS and SAT-1B programs is that the MAIN routine in SAT-IB also performs 
the operations of the subroutine BOZO (in SAT-lS). The GRAPH routines are 
essentially identical between the programs. 
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SECTION 4 


SATELLITE- IV PROGRAM 


The discussions in SECTION I for items 1.1, 1.2 are equally valid for all 
SAT-l( ) programs. 

EXAMPLES, FLOW CHARTS, LISTING: One sample test problem is executed by 

the SAT-1V program. Since most of the errors have been previously 
discussed, the current problem exemplifies a successful submission. 
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The listing of the SAT-1V program follows. The flow charts for this program 
are the same as for the SAT-1B program ^discussed previously. 
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CALL CHSIZV 12 , 2 ) 

CALL RITSTV ( 13 , 19 , TABL l V ( 


IF (ZZ< JI.LE.YMIN) YMIN = ZZ(J) 

IF (ZZ(J).GT.YHAX) YMAX = ZZ(J) 

IF IRRAOI JKGT.XHAX) XHAX = RRAD(J) 
400 CONTINUE 
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